Non local screening in a vortex line liquid 
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We show that the recent experiments Q reporting the onset of the non-local conductivity in the 
vortex state of Y BaCuO single crystals indicate the presence of a new liquid phase of vortices. This 
phase is intermediate between the normal metal and the Abrikosov lattice. We use the mapping of 
the vortex problem to the problem of bose liquid to determine theoretically the properties of the 
proposed vortex liquid phase and compare them with the data. 
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In a recent Letter H. Safar et al reported 0] the re- 
sults of the transport measurements in the mixed state 
of YBaCuO single crystals with the field applied along 
c-direction of the crystal performed in a transformer ge- 
ometry shown in Fig 1. Two sets of measurements were 
done. In the first one the current was injected through 
the pair of contacts (1,4) and the potential was mea- 
sured between contacts (2,3) or (6,7). In the second set 
the current was injected through contacts (1,5) and the 
potential was measured between contacts (2, 6), (3, 7) or 
(4, 8). In the first setup the current flows predominantly 
along ab plane of the crystal whereas in the second setup 
the current flows mostly in c direction. 



Fig. 1. Geometry of dc transformer experiment. 
Heavy dots denote contacts. 

It was found JO that in the first setup the potential 
drop between contacts (6, 7) is smaller than the drop be- 
tween (2, 3) at high temperatures, but becomes undistin- 
guishablc from it below some critical temperature T t h- 
If interpreted in terms of the apparent resistivity ratio 
Pc/ Pah such observation implies that this ratio tends to 
zero at T — > T t h(H). This temperature T t h(H) is sig- 
nificantly higher than the transition temperature T g (H) 
defined by p a b(T g ) = 0. However, the experiments of 
the second setup showed that the ratio of the potentials 



V48/V26 becomes much larger when T — » T t h(H) indicat- 
ing that apparent ratio p c / p a b extracted from these mea- 
surements tends to infinity. Finally, an independent set 
of measurements indicate that resistivity in c-direction 
extracted from measurements with the uniform current 
p c (T th )=0 §. 

Safar at al interpreted their results as an evidence for 
non local transport in the mixed state; in this note we 
shall show that this non local transport may be a sig- 
nature of a phase transition into the novel vortex liquid 
state characterized by a non-zero phase rigidity in the 
direction of the magnetic field HQ . We shall argue that 
T t h should be identified with the transition temperature 
into the intermediate vortex liquid state. 

We begin with the estimates of the parameters of the 
vortex system in the regime of the experiments jj],^). 
We use the value of the penetration length A = 1400A 
and anisotropy ratio M/m = 50 to estimate a vor- 
tex entanglement length (cf. ||) in c-direction: L = 
1.5(1 — t)/B pm where t = T/T c and B is measured 
in Tesla which gives L ~ 0.03 pm in the conditions of 
HH. Thus, the vortices are entangled on a much shorter 
scales than the sample thickness (0.03 mm). This esti- 
mate rules out the simple explanation of the experiment 
|l| that each individual vortex remains straight on the 
scales of the sample thickness. The absence of current 
dissipation in c-direction shows that although vortices 
wander at short scales they remain relatively straight at 
large scales due to collective effects. The dissipation of 
the uniform in-plane current shows that the vortex lattice 
is not formed. 

In order to describe the liquid of vortices state quan- 
titatively we mapped the vortex problem to the model 
of strongly interacting two dimensional bosons We 
shown that the problem of Gibbs equilibrium state of the 
system of infinitely long vortices at temperature T is for- 
mally equivalent to the quantum ground state problem 
of two-dimensional interacting Bose liquid (with the " 2D 
Plank constant" T12D = T) at zero temperature. The 
direction of the vortex space along the magnetic field 
(z || c) plays the role of imaginary time of the 2D quan- 
tum theory. The interaction between bosons in the quan- 
tum model has two parts: static 2D Coulomb part and 
dynamic transverse current-current interaction mediated 
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by the fluctuating 2D "magnetic" field (cf. §,|). This 
mapping identifies the Abrikosov lattice and the crys- 
tal of bosons, the normal metal and the Bose superfluid 
liquid. In we also found a new vortex liquid phase 
between Abrikosov lattice and normal metal which is 
mapped into the "normal" (i.e. nonsuperfluid) liquid of 
bosons. 

Here we shall use this mapping to find the response 
of the vortex phase in the conditions of experiment |Q . 
In we derived the duality relations between the re- 
sponse functions in the vortex system and the response 
function for the boson system and found that the new 
intermediate phase has non-zero superfluid density in c- 
direction: j c = —nfA c but has a finite resistivity to 
a uniform current in a6-directions. The phase transition 
line normal metal-new phase predicted in Ref. |^|, Fig. la 
is surprisingly close to T t h(H) line reported in [^]J^]. Be- 
low we show that in this state the non-uniform in-plane 
current which satisfies J jabdz — is non-dissipative. 
This resolves the apparent contradiction of observations 
because in the set up yielding p a b/p c — > the in-plane 
current obeys J j a bdz = 0. 

We prove the existence of a non-dissipative current 
jabilz 7^ 0) repeating the derivation of response func- 
tions (cf. Section VI. A. 1 of Q) for the in-plane currents. 
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where n || (q) is the longitudinal response of the dual bose 



liquid and g 2 = 



4ttA 2 



is its interaction constant i.e. the 



analog of the electric charge e 2 D ■ Here and below we ex- 
plicitly consider only isotropic superconductors; we shall 
restore anisotropy factor m/M only at the very end using 
general scaling arguments ||. 

The intermediate phase corresponds to the normal liq- 
uid of bosons. Since q z play the role of the frequency in 
the bose model Xln (q) is strongly dependent on the ratio 
q z /q±. In the limit q± — > the longitudinal response of 
the boson system is described by the finite effective 'con- 
ductivity' a: (/ 2 n||(<7 2 ) = u\q z \. To estimate the effective 
conductivity we note that normal bose liquid is realized 
in the regime of the strong interaction when the dimen- 
sionless interaction parameters are the order of unity jjj . 
In this regime the only combination of parameters with 
dimensionality of conductivity is the "quantum" conduc- 
tivity ctq = e 2 D /h 2D . Identifying g — > e 2D , T — > h 2Dl we 
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estimate a ~ gL. Using this expression for H||(g 2 ) we 
get the correlator of the in-plane electromagnetic vector 
potential: 



D(q z ) = 



4ttT 
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Eq. (g) shows that the in-plane magnetic field decays 
at large distances (z ^> A) as B(z) ~ B A in this state, 
being screened on a scale of A by a non-dissipative cur- 



rent j 



along the boundary. To see this, con- 



sider the effect of the external current flowing in the x 
direction on the edge of the sample j x xt (z) — J ext 5(z). 
According to Eq.(|2|) it induces magnetic field B y (q z ) = 
q z ■ D(q z ) ■ J ext oc l/(q z + sgn{q z )A~ 1 ). After integra- 
tion over q z one gets the induced field B y (z) ~ Bqj and 
the shielding current j x {z) oc dB y /dz. Note that in the 
superfluid phase of bose-liquid one would get a oc l/q z 
and H||(g 2 ) ~ const, so the electromagnetic correlation 
function would acquire its usual for normal metal form 
D(q z ) oc q z 2 . The universal length scale A is about 
400 fim at T — 90 K which is larger than the sample 
thickness d ( note that the rescaling taking into ac- 
count anisotropy does not affect the relation between A 
and d, since both quantities scale in the same way) . The 
existence of non-dissipative in-plane current at q z ^= 
in the vortex liquid implies that p xx (q z =/= 0) = at 
T < T t h(H). We interpret observations as a sig- 
nature of vanishing of p X x(qz 7^ 0) at T — > T t h{H) + , 
which supports a phenomenological description proposed 
in Ref. §. 

The Eq.fQ) for D(q z ) was derived in the limit of in- 
plane homogeneous current, i.e. q± — * 0. Using the 
conventional diffusion form for the density-density corre- 
lator 
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where A = X 2 a = 



we see that Eq.(||) remains valid if q\ <C q z /T^- Here T> 
has the meaning of "diffusion coefficient" associated with 

"conductivity" a. Thus, a w e% D V C^f-J where the 
last factor in parenthesis is the density of states for 2D 
particles with mass m 2 D — (<I>o/47rA) 2 . The "quantum- 
limit" expression for cr, which we used above means that 
the mean free path of bosons is of the order of their sep- 
aration, i. e. T> ~ n 2 r>/m 2 r)- 

Now we estimate the relative size of wavevectors q z 
and q\ for the conditions of the experiment M . After the 
rescaling || takes into account with the mass anisotropy 
factor m/M = 0.02, the relevant wavevectors become 
q z w yJm/M j d, q± « l/d p i ane , with the sample thick- 
ness d — (3 — 6) • 10~ 3 cm and the relevant lateral dimen- 
sion dpiane ~ 0.05cm. With the above estimates for T>, 
the condition q\ -C q z /T> is satisfied for any reasonable 
sample size. 

We expect also that the results of the dc transformer 
measurements in thick samples (d > A) should be quali- 
tatively different from the samples reported in QJ^] due 
to the screening effects. Finally, we note that the re- 
sults obtained in twinned samples for the tilted field will 
be severely modified by the interaction between vortices 
and twins because in these conditions each vortex line 
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intersects the twin boundary which affects strongly the 
motion and equilibrium positions of vortex lines. 

In conclusion, the observations jlj,0 strongly support 
the existence of a new vortex line liquid state sandwiched 
between the normal state and the vortex glass. 

We are grateful to H. Safar for useful discussions. We 
acknowledge the support of the innovation partnership 
grant of the State of NJ and of the grant M6M000 from 
International Science Foundation. 



[1] H. Safar, P. L. Gammel, D. A. Huse, S. N. Majumdar, L. 

F. Schneemeyer, D. J. Bishop, D. Lopez, G. Nieva and F. 

de la Cruz, Phys. Rev. Lett. 72, 1272 (1994). 
[2] F. de la Cruz, D. Lopez and G. Nieva (unpublished). 
[3] M. V. Feigelman, V. B. Geshkenbein and V. M. Vinokur, 

JETP Letters, 52, 546 (1990). 
[4] M. V. Feigelman, V. B. Geshkenbein, L. B. Ioffe and A. I. 

Larkin, Phys. Rev. B 48, 16641 (1993). 
[5] D. R. Nelson Phys. Rev. Lett. 60, 1973 (1988). 
[6] G. Blatter, V. B. Geshkenbein, A. I. Larkin, Phys. Rev. 

Lett. 68, 875 (1992); cf. also detailed discussion in G. 

Blatter et al, Rev. Mod. Phys. (1995). 
[7] D. A. Huse and S. N. Majumdar, Phys. Rev. Lett. 71, 

2473 (1993). 



3 



